Dynamical evolution of correlated spontaneous emission of a single photon from a 

uniformly excited cloud of N atoms 
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We study the correlated spontaneous emission from a dense spherical cloud of A'^ atoms uniformly 
excited by absorption of a single photon. We find that the decay of such a state depends on the 
relation between an effective Rabi frequency Q oc \/iV and the time of photon flight through the 
cloud R/c. If QR/c < 1 the state exponentially decays with rate Q^R/c and the state life time is 
greater then R/c. In the opposite limit Q,R/c ^ 1, the coupled atom-radiation system oscillates 
between the collective Dicke state (with no-photons) and the atomic ground state (with one photon) 
with frequency il while decaying at a rate c/ R. 



The problem of a single photon absorbed by a cloud 
of N atoms followed by correlated spontaneous emission 
is a problem of long standing interest. Dicke [1] first 
noted that the radiation rate from a small dense cloud 
is abnormal. In his words "... the greatest radiation 
intensity anomaly occurs in the transition to the ground 
state" [2]. In particular he showed that the collective 
decay rate for the symmetric state with one excitation is 
r^v = N^ (as indicated in Fig. Ilia)), where he assumed 
the atomic volume to have dimensions small compared 
with the radiation wavelength. 

In our present work, we are interested in the time evo- 
lution of a specially prepared state obtained by absorp- 



tion of a single photon [3J|31[S]. We report novel dynam- 
ical oscillations in the evolution of the quantum state of 
the atom cloud, even without the existence of a cavity. 
Fig. [T]summarizes the main results of this paper. It is as 
if the atomic cloud acts to form a new "cavity" with the 
atom cloud volume V replacing the virtual photon vol- 
ume Vph defined by the electromagnetic cavity; that is the 
usual vacuum Rabi frequency fivac = (p/^) '\/^w/(eoVph) 
is replaced by fig = {p / h) y^hw/JeoV) in the present 
problem, where p is the electric-dipole transition ma- 
trix element, hw is the photon energy and eo is free space 
permittivity. 
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FIG. 1: Comparison of the different dynamical behavior for the correlated spontaneous emission from an A atomic cloud 
described by the |+)ko state of Fig. 2b. The single atom spontaneous decay life time is taken to be tq = 10 ns (7 = l/2ro). We 
assume that atomic density is lO^^cm"^ and the resonant photon wavelength is A = 1/xm. Plot (a) corresponds to the case when 
cloud radius is equal to A/2, hence the number of atoms is Na = 5240, then the state decay time is Ta ~ ro/Na ~ 1.9 x 10~^ 



ns. In plot (b) the cloud radius is r = lOA, Ni, = 4.2 x 10^, 
the atomic cloud is i? = 1 mm which yields Ac = 4.1 x 10^ 
27r/n = 0.74 X 10"=^ ns. 



Tb = Z2-K^r^To/27Nty = 2.7 x 10""* ns. In plot (c) the radius of 
8R/6c = 4.4 X 10~^ ns, while the period of oscillations is 



Here we analytically solve the equation of motion for 
the system in two regimes. In one regime, we find that, 
for finite atom cloud size, the atomic excitation will de- 
cay with a rate determined by the photon escape time, 
together with fast dynamical oscillations (as indicated in 
Fig. [iTc)) with the effective Rabi oscillation frequency fi 
proportional to \/N . 

However, if the oscillation period is much greater then 
the time of photon flight through the cloud, the atomic 
state will decay exponentially. For a small cloud {R ^ A) 
the decay rate is N^, where 7 is the single atom de- 
cay rate, see Fig. la. For a larger cloud (i? » A but 
^R/c <C 1) the decay rate goes as Nj{X/R)'^ as in Fig. 
lb. Finally, for an even larger cloud {Q,R/c ^ 1) the 
probability of photon emission oscillates while decaying 
with a rate ~ c/ R, see Fig. [it. 

These two regimes (b and cj are determined by whether 
the system persists memory effect, i.e., non-Mar kovian or 
Markovian regimes. We discuss connection with previous 
work [1 [71 [HI H [M EI] at the end of this paper. 

We consider a system of N two-level (a excited and b 
ground) atoms, initially one of them is in the excited state 
a (with no information which one), Ea — Ej, ^ huj, and 
the multi-mode radiation field is in the vacuum. Atoms 
are located at positions Vj {j = 1, ...,N). The whole set 
of states can be expressed as those in Fig. 2 |T, where 



\j) = \bi,b2, ■■■,bj-i,aj,bj+i...bN) represents the state 
in which the j*^ atom is exited but the others are in the 
ground state and \g) — |6i, 62, ..., b[^) is the state with all 
the atoms in the ground state. The atomic state prepared 
by uniformly absorbing one single photon with wavevec- 
tor ko is exactly the |+)ko state of Fig. 2b. In the limit 
-R S> A we focus here the |-f )ko state is approximately an 
eigenstate of the system [4j . This makes the basis set of 
Fig. 2b preferable. The state vector for the atom-field 
system at time t can be then written as 

\nt)) = [/3+(<)l+)ko + /3i(t)|l)ko + /32(t)|2)ko + - 



-K/3^_i(i)|7V-l)ko]|0) + 5]7k(t)|5)|lk), (1) 



with initial conditions /3+(0) = 1 and all other probability 
amplitudes are zero. In the dipole approximation the 
atom-field interaction is described by the Hamiltonian 



N 

Hint = X! X! ^9k cTjal exp{i{u!k - ^)t - ik • r^) + adj 
k j=i 

(2) 
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FIG. 2: Timed Dicke states associated with absorption of radiation of wave vector ko: (a) the initial state |+)ko decays directly 
to the graund state \g). (b) the timed Dicke states corresponding to single photon excitations. 
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where aj is the lowering operator for atom j, Sk is the 
photon operator and g^ = p\/ fu-Ok / {f-oVph) is the atom- 
photon coupling constant for the k mode, u>k is the pho- 
ton frequency and lo = ck^ is the energy difference be- 
tween level a and b, c is the speed of light. For simplic- 
ity, we neglect the effects of photon polarization. The 



dynamical evolution is then totally determined by the 
Schrodinger's equation. Let us consider first the two 
atoms problem, and call the state |l)ko = l^)ko- The 
state vector for the two atoms plus field system is then 



given by 

I^W>-/3+WI+)ko|0)+/3-(<)|-)kol 



-E^k(t)|5)|lk) 



(3) 
As is shown in Rcf. [T, the probability amplitude /?+ 
and /3_ are coupled due to the fact that they decay to a 
common ground state, that is 



/?+ = -7++/5+-7+-/3-, 



/?_ = -7-/3- - 7-+/34 



(4) 



(5) 



This was first pointed out by Agarwal jT2]. It is closely 
related to earlier work of Fano, and is often referred to 
as Fano coupling. We call it Agarwal-Fano coupling. 

However, when we consider a sphere with N atoms, 
where iV 3> 1, and R ^ X there is no Agarwal-Fano 
coupling. Instead we now find ^ 

1 /•* ^ 



k JJ=1 



.exp[z(k-ko)-(r,-r,)]/3+(t'). 



(6) 



For a dense cloud one can treat the atom distribution 
as continuous, we then have J2i ~^ i^/^) J d-'r f dr' , 
where V = 47ri?^/3 is the volume of the spherical atomic 
cloud. The summation over k can also be replaced by in- 
tegration ^j^ — + Vp/i/(27r)'^ / dk, where Vph is the pho- 
ton volume. Then the equation of motion reads 



where fi = v -^f^o is an effective Rabi frequency. There- 
fore in the limit i? —s- oo the atomic state undergoes 
harmonic oscillations with the effective Rabi frequency 



p+{t) ^cos{nt). 



(10) 



To find a solution of Eq. (l7| at finite R, but yet koR ^ 
1, we rewrite it as 

/3+(t) = -?%^ [dkgl ^c^i'/340e^^---)(*'-*^5(k,i?)^ 



nV^ 



where 



(11) 



S{k,R) = — / drexp[i(k-ko)r] 
47r Jy 



sin(|k-ko|i?) i?cos(|k-ko|i?) 



|k-ko|3 



|k-ko| 



(12) 



Next we approximate g^ ~ g\ and replace integration 
over k by integration over p = k — kp. The main con- 
tribution to the integral comes from the region p < \/ R. 
That is under the exponent one can replace k — k^) "^ 
ko • p/fco. Then Eq. (Ill reads 



P+it) 



N^l I dp / dt'f3+{t')e 



''ipb'cko-p(t'-t)/fco] 



sin(pi?) Rcos{pR) 
Integration over directions of p yields 



(13) 



/3+(i) = - 



Vph N 
(27r)3 V^ 



dk / dr 



dv'gl / dt'p+{t') 





P+{t) = —Nnl / pdp dt'f3+{t') 
c Jq Jo 



■ exp[i{ujk - oj){t' ~t)+ i{k - ko) • (r - r')] . (7) 

In the limit i? — » cx) integration over r' gives the delta- 
function 

/"dr'exp[-i(k-ko)r'] = (27r)=^(5(k - ko), 



and thus we obtE 



p+{t) = -Nnl / dt'(3+{t') 

Jq 



(8) 



where we have defined fio — {p / h) ^ hw / {eoV) , which is 
like vacuum Rabi frequency but with atomic volume V 
replacing photon volume Vph- 

Differentiating both sides of Eq. (l8| yields a harmonic 
oscillator equation 



^+(t)+f}2^+(i) = 0. 



(9) 



sin[cp(i' - t)] 
{f - t) 



sin {pR) i? cos {pR) 



po 



pz 



To integrate over p we use the following formula 



(14) 



sin[cp(t' - i)] 
'"' (t'-t) 



sin (pi?) R cos {pR) 



■KC 

96 



[I6i?3 + 12ci?2(i' - t) - c^{t' - tf] , c\t' -t\<2R 
0, otherwise 

(15) 
which gives 

/3+(t) = -^Nnl J dt'l3+{t') [l6 + 12 -^(i' - t) 



I .^3 



^3 



(t'^t) 



e [c(i' -t)+ 2R] 



(16) 



Next we note that the function 16+12 j^{t'-t)-^{t'-tf 
and its derivative over t' is equal to zero when c(i' — t) + 
2R — 0. Taking derivative of both sides of Eq. (IT6| twice 
we obtain: 



P + it) 



-Nfll 



P+it) 



3c 



P+it) 



3c' 



.3 pt 



SR^Jo 
Next we assume that 



dt'P{t'){t' - <)e [c{t' -t) + 2R] 



nR 



Nft'o^ > 1, or 

c^ c 



>1. 



(17) 



(18) 



Then one can omit the last term in Eq. ( 17) which yields 



3rn2 

(3^{t) + n'(3+it)-^^P+{t)^0. (19) 



Solution of Eq. ( 19 ) under the condition ( 18 1 is given by 

3c 



/3+(t) — cos(rii) exp 



8R 



tU 



(20) 



which describes rapid oscillations with the effective Rabi 
frequency fl superimposed by the exponential decay. The 
state decays during the time of the photon flight through 
the atomic cloud. The emitted photon is reabsorbed and 
reemitted many times before it leaves the cloud. 

In the opposite limit, flR/c <^ 1, one can use the 
Markovian approximation. We integrate Eq. (16) over t' 



assuming (3+{t') is a slow varying function of t' and ap- 
proximate /3+(i') ~ P+{t)- Then for t > 2R/c we obtain 



$+{t) = -r/3+(i), 

which yields an exponentially decaying solution 



/?+(t)=/3+(0)e 



-rt 



(21) 



(22) 



Here T = 3n^R/4c = 27N-f/8{kaR)^ and 7 = 
{uj^p'^)/{6TTeohc^) is the spontaneous decay rate for one 
atom. 

Discussion: Similar problems have been investigated 
in the past several decades. Cummings fP. considered 
the spontaneous emission of a single atom which is ini- 
tially excited in the presence of iV — 1 initially unexcited 
identical atoms, when there are M accessible radiation 



modes. He showed that such an extended system oscil- 
lates between the ground sate and the excited state with 
an effective Rabi frequency fi ~ y/N. Such modification 
of the spontaneous emission of one atom in the presence 
of TV — 1 atoms inside a cavity has been studied since 
then [8], [9]. Buzek [10] studied the dynamics of an ex- 
cited atom in the presence of TV — 1 unexcited atoms 
in the free space and predicted that there is a radiation 
suppression but did not report dynamical oscillations. 

The effective Rabi frequency fi = VND,o we found 
from quantum mechanical consideration can be written 
as r^ = v^(3/47r2)7cj(iV/F)A3 := y^mJ^Pje^. This re- 
sult is analogous to the plasma frequency and can be ob- 
tained in a classical model by treating atoms as classical 
harmonic oscillators [T5^. Indeed, replacing the electric- 
dipole transition matrix element by p = e ■ d, where 
d — ^Jhjnvxi is the oscillator length, yields precisely the 
plasma frequency fi = ^ne^ jmt^. 

Relevant experiments have been carried out by the 
groups of Lukin [14] , Kuzmich [15] , Kimble [16] , Vuletic 
[T7], Harris [TB] et al. [TH]. For realistic physical situa- 
tions such as n = lO^^'cm-^, w/27r = 6 x lO^^^Hz, i? = 10 
cm (iV = 4 X lO^'^) and \p\ = lO'^^C-m, we obtain that 
the state decay is accompanied by a few oscillations with 
the effective Rabi frequency O « 2 x 10^°Hz and the de- 
cay time is about R/c ^ 3 x 10~^°s. One can observe a 
crossover to the exponentially decaying regime, e.g., by 
decreasing the size of the atomic cloud. 

In summary, we study correlated spontaneous emission 
of a totally symmetric A^-atom state prepared by an ab- 
sorption of a single photon. This is an extension of the 
result obtained in Refs. [31 H]. Decay of such a state oc- 
curs via photon emission in the direction of the incident 
photon for large enough density. We found that time evo- 
lution of the initial state depends on the relation between 
an effective Rabi frequency ft ex \/Ngkg and the time of 
photon flight through the cloud R/c. If ilR/c < 1 the 
state exponentially decays with the rate il^R/c which 
is determined by the Dicke superradiance rate Nj re- 
duced by the factor of \^ /R? due to smaller finite state 
phase volume in the case of directional emission. In the 
opposite limit ^R/c :$> 1 the decay is accompanied by 
oscillations with the effective Rabi frequency fi and the 
decay time is given by R/c. 
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